Vibrational edge modes in intrinsically inhomogeneous doped transition metal oxides 
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By applying an unrestriced Hartree-Fock approximation and a Random Phase approximation to 
multiband Peierls-Hubbard Hamiltonians, we determine the phonon mode structure in models of 
transition metal oxides in the presence of intrinsic nanoscale inhomogeneities induced by hole doping. 
We identify low frequency local vibrational modes pinned to the interface between regions of distinct 
electronic structure and separated in frequency from the band of extended phonons. A major 
characteristic of these "edge" modes is that their energy is essentially insensitive to the doping level, 
while their intensity increases with the density of interfaces (and thus doping level). We argue that 
the presence of such modes is a typical feature of systems with phase separation, including cuprates, 
nickelates, manganites, and bismut hates. We specifically address the experimental signatures of 
these modes in lattice inelastic neutron scattering. 



Recent advances in experimental techniques provide 
strong indications that a broad class of electronic ma- 
terials, commonly referred to as strongly correlated, ex- 
hibit intrinsically inhomogeneous electronic phases [1-7]. 
Strong interactions in these materials drive phase sepa- 
ration between regions of distinct electronic structure [8] . 
However, a competing interaction or electronic kinetic en- 
ergy can lead to frustration of the global phase separation 
and to formation of nanoscale inhomogeneities, which can 
take the form of, e.g., stripes [9] or clumps [10]. It has 
been realized early on that unlike the conventional spin 
and charge density wave phases, the ordering of intrin- 
sic inhomogeneities is sensitive to thermal and quantum 
fluctuations [11], as well as static disorder. This makes 
an unambiguous identification of such phases, and their 
distinction from the more common density wave instabil- 
ities problematic, and therefore requires development of 
new diagnostics, which are more sensitive to the unique 
structure of such phases. In this Letter, we propose one 
such diagnostic tool which relies on detection of localized 
edge modes pinned to the interface between the regions of 
different electronic structure. While the edge (interface) 
modes are predicted to appear in every degree of freedom 
coupled to the inhomogeneity (charge, spin, lattice) [12], 
here we focus of the lattice vibrational edge modes which 
lend themselves readily to study with inelastic neutron 
scattering. 

We consider two particular classes of transition metal 
oxides that have attracted significant attention recently 
- cuprates and nickelates. The former is the material 
exhibiting high-T c superconductivity, and the latter is a 
related compound with Cu replaced by Ni. The Ni-based 
compounds are typically insulating, with unambiguously 
demonstrated stripe order [1]. In cuprates, on the other 
hand, the presence of stripes and their possible influence 
on superconductivity are subjects of intensive debate. 
Since it appears that stripes can significantly affect super- 
conductivity [13-15], it is important to establish whether 



they are indeed present in the material. Our purpose 
here is to elucidate typical manifestations of stripe order 
through a comparative theoretical study of cuprates and 
nickelates. In particular, we compare these theoretical re- 
sults to observations of anomalous phonon modes in the 
cuprates [16-18] and nickelates [19,20]. The results also 
apply to edge mode signatures in other doped transition 
metal oxides (and related complex electronic materials), 
as found in recent experimental studies [19,21-23]; equiv- 
alent modes have been predicted and identified in conju- 
gated polymers and quasi- ID (one-dimensional) charge 
transfer solids [24,25]. 

We first analyze the signatures of the edge modes in 
neutron scattering for doped nickelates. To model the 
MO2 planes, we use a 2D four band extended Peierls- 
Hubbard Hamiltonian, which includes both electron- 
electron and electron-phonon interactions [26,27]: 
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Here, c\ m(T creates a hole with spin a and site i in orbital 
m (there are two nickel orbitals (Ni d x 2_ y 2 and d 3z 2_ 1 ) 
and two oxygen orbitals (O p^ and p y ) included). The 
bare Ni-0 hopping t irn j n has two values: t p d between 
d x 2_ y 2 and p and ±t p d/V3 between d^ z 2_i and p. The 
O site electronic energy is e p and the Ni site energies 
are and + E z , with E z the crystal-field splitting on 
the Ni site. The correlation part of the Hamiltonian is 
modelled by 

He' = ( U + 2 J ) n imrnimi ~ Yl 2JS im .S in (2) 
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Here, U is the on-site Ni Coulomb repulsion, J is the 
Hund coupling that favors a high spin state on Ni sites, 

/, with a the vector of 
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Pauli matrices. The electron-lattice interaction causes 
modification of the Ni-0 hopping strength through the 
oxygen displacement mf t irn j n (uij) = t irn j n (l ± 
where +(— ) applies if the Ni-0 bond shrinks (stretches) 
for a positive it also affects the Ni on-site energies 
£m(uij) = £d + P^2j (=^%')> where the sum runs over 
the four neighboring O ions. The other oxygen modes 
couple to electron charge more weakly, and are neglected 
here for simplicity. In our calculations we use the follow- 
ing set of parameters [28,29]: e p — — 9, U = 4, J = 
1, E z = 1, and K = 32/ A (all in units of t v d\ we take 
t p d = 1 eV). The electron-lattice couplings a = 3/ A and 
(3 = 1/A were varied to best fit the neutron scattering 
data. To approximately solve the model, we use unre- 
stricted Hartree-Fock (HF) combined with an inhomo- 
geneous Random Phase Approximation (RPA) for linear 
lattice fluctuations [26,29] in a supercell of size N x x N y 
with periodic boundary conditions. In this model, doped 
holes tend to segregate into the stripes due to compe- 
tition between magneto-elastic interaction that favors 
global electronic phase separation and electronic kinetic 
energy that favors uniform carrier density. It is worth 
noting that other competing interactions can produce 
stripes, clumps and other inhomogeneities [30,31], but 
the local coupled charge-spin-lattice dynamics governing 
edge modes on these templates can still be modelled with 
the above Hamiltonian and RPA analysis. 

The output of the calculation is the inhomogeneous HF 
groundstate and the phonon eigenfrequencies and eigen- 
vectors. From the phonon eigenmodes, we calculate the 
corresponding neutron scattering cross section, defined 
as 
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where R/(t) = + di + u/(t) is the position of the 
l-th oxygen atom expressed in terms of the location of 
the unit cell origin R^, position within the unit cell d/, 
and time-dependent vibrational component u/(t). For 
phonon modes with u/(t) oriented along the correspond- 
ing metal-oxygen bonds, on the O x sublattice di = |f 
and ui = x\x , and on the sublattice di = and 
= yiy. The scalar displacements can now be expressed 
in terms of the normal modes z n as xi(t) = ^ n a xi,nZ n (t) 
and yi(t) = V &yunZn(t) . Making first order expansion 
in the oxygen displacements, we obtain 



^(k,c) = ^feK, n | 2 + ^K, 
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the thermal population of the phonon mode n. This is 
a generalization of the usual neutron scattering intensity 
expression [32] for the case of phonons with a larger real 
space unit cell. In this paper we plot *S(k, o;)/|k| 2 for k- 
directions sampling longitudinal modes, consistent with 
the common experimental convention. 
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FIG. 1. Neutron scattering spectra and the ground state 
structures for the undoped (left panels), and hole-doped 
(x = 0.2) MO2 planes. The neutron scattering spectra are 
symmetrized with respect to directions along and across the 
stripe; in all plots q runs from to 2ty. The relative intensity 
color scheme is defined in Fig. 3. In the lower plots, black 
circles represent Ni sites and red circles represent O. The ra- 
dius of a circle is proportional to the corresponding site hole 
density. Arrows centered on circles show the magnitude and 
direction of spin. The green lines originating from O sites 
indicate the magnitude and direction of the equilibrium O 
displacements due to the formation of a stripe. 

The ground state configurations and the corresponding 
neutron scattering intensities for MO2 planes are shown 
in Fig. 1. The ground state configuration in the un- 
doped case is a simple antiferromagnet (AF), while, in 
the doped (x = 0.2) case, the holes spatially segregate 
into a Ni centered stripe along the diagonal of the NiC>2 
plane. The neutron scattering intensities are given in the 
original Brillouin zone, symmetrized with respect to the 
orientation of the stripe. 
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The most striking new feature in the dynamic struc- 
ture factor that appears upon doping is the creation of 
the low frequency modes (u ~ 75 meV) that split off from 
the main branch. These low frequency modes have been 
observed experimentally with inelastic neutron scatter- 
ing [19,20]. We find that with increasing doping (e.g., 
x = 0.5) the spectral weight in the new branch increases. 
However, its energy does not significantly change, nor is 
there a dependence of the spectrum on the stripe spac- 
ing. Intensity to create the local modes is taken from 
the extended phonon mode branch, with accompanying 
phase shifts of the extended phonons. There is no simple 
dispersion that can be associated with the local modes 
and their wave vectors are related to the characteristic 
spatial extent of the eigenmodes. This is consistent with 
the doping insensitivity of the phonon "anomaly" ob- 
served in neutron scattering experiments for nickelates 
reported in Ref. [20]. The presence of the phonon edge 
mode with doping-independent energy is similarly con- 
sistent with the doping-insensitivity of the "kink" energy 
in Angle Resolved Photoemission Spectra (ARPES) re- 
ported recently [33] for cuprates (discussed later). 

In order to elucidate the structure of these new modes, 
in Fig. 2 (left panel) we plot a representative mode from 
this branch. Clearly, this is a mode localized around the 
stripe edge. For comparison, in the right panel we plot 
a mode from the high energy branch. Although modi- 
fied (phase-shifted) due to the presence of a stripe, this 
mode preserves its original spatially extended structure, 
characteristic of the undoped case. 
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FIG. 2. Representative vibrational eigenmodes for the case 
of a x = 0.2 doped N1O2 plane. Black and red circles represent 
the equilibrium positions of Ni and O atoms, respectively, 
and green lines show the magnitude and direction of the O 
displacement in the mode. On the left, a low frequency mode 
localized in the vicinity of the stripe (E = 75.1 meV); on the 
right, an extended phonon with energy E = 83.4 meV. 

We now turn to the case of doped Cu02 planes for com- 
parison. To this end, we analyze the three-band model 
analogous to Eqs. (1) and (2), but with only one orbital 
on the metal site (Cu d x 2_ y 2). The electron-electron in- 
teractions that we include are the Cu on-site repulsion 
Ud = 8, the O on-site repulsion U p = 3, and the nearest- 
neighbor Cu-0 repulsion U p d = 1. The remainder of the 
parameters are e p — €d — 4, K = 32/ A, a = 4. 5/ A and 
/3 — 1/A, with t p d = 1.2 eV [12]. Results of calculations 
for undoped and doped cuprates are presented in Fig. 3. 
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Neutron scattering spectra and the ground state 
for the undoped (left panels), and hole-doped 
(x = 0.2) CuC>2 planes. The neutron scattering spectra are 
symmetrized with respect to directions along and across the 
stripe. In the lower plots, black circles represent Cu sites and 
red circles represent O. The radius of the circle is proportional 
to the corresponding site hole density. Arrows centered on cir- 
cles show the magnitude and direction of spin. The green lines 
originating from O sites indicate the magnitude and direction 
of the equilibrium O displacements due to the formation of a 
stripe. 

While generally similar to the MO2, in the Cu02 case 
doping produces two low-frequency branches of localized 
modes that split off from the band of extended phonon 
modes. Representative modes from each localized branch 
are shown in Fig. 4. The lower frequency band centered 
around ~ 40 meV (left) corresponds to the oxygen dis- 
placements along the stripe and the high frequency one at 
~ 73 meV to the vibrations perpendicular to the stripe. 
(The case shown in Fig. 4, right, is a localized asym- 
metric breathing mode [16,17].) Experimental evidence 
shows that mode softening of about 15 percent occurs 
with doping in the cuprates, consistent with the soft- 
ening predicted here for the perpendicular edge mode 
[16,17]. Recent neutron measurements in twinned and 
detwinned YBa2Cu306+5 [18] have attempted to resolve 
the phonon spectrum both parallel and perpendicular to 
the purported stripe direction [34] , but conclusive identi- 
fication has been hampered by the degeneracy and mix- 
ing of many modes near the zone boundary. As yet, there 
has been no search for the parallel edge modes which are 
predicted here to have a huge softening of ~ 60 percent. 
However, there are distinct anomalies in the longitudi- 
nal phonon spectrum of I^-^Sr^CuC^ at 30 meV that 
have not been studied in detail [17]. Also, it is interesting 
to mention a possible connection to the "magnetic res- 
onance" mode reported at 41 meV [35]. Note that this 
low-energy mode (see Fig. 4, left) can only exist for the 
vertical stripe and hence is absent for the diagonal stripes 
in the case of MO2 (Fig. 1). 
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FIG. 4. Localized vibrational eigenmodes for the case of a 
x = 0.2 doped C11O2 plane. The stripe is centered along the 
middle vertical row of Cu. Unlike the case of diagonal stripes 
in Ni02, there are two localized branches, one that corre- 
sponds to the oxygen vibration parallel to the stripe (low fre- 
quency, E = 31.2 meV, on the left), and one that corresponds 
to the oxygen displacements perpendicular to the stripe (high 
frequency, E = 73.3 meV, on the right). 

Despite the impressive agreement between the energy 
scales of the predicted and observed modes, it is impor- 
tant to note some deviations from the reported exper- 
imental data in terms of the mode intensities and dis- 
persions. We speculate that these discrepancies may be 
due to the detailed local oxygen environments: either (a) 
more complex intra-unit cell distortions, i.e. local oxy- 
gen polarizations (e.g., from local confinement or pair- 
ing effects), which would require more detailed modelling 
than the constrained O displacements used here, or (b) 
more complex stripe configurations than the regular lin- 
ear stripes we have considered above. (Such textured 
stripe patterns are indeed found in certain models for 
the origin of stripes [8,31]). 

We conclude that the formation of split-off local 
modes, well known in the polaronic physics of conjugated 
polymers and quasi-lD charge transfer salts [36,24,25], is 
also expected to be prominently manifested in the case 
of other doped strongly correlated materials, including 
nickelates, cuprates, and bismuthates [1-7]. Elsewhere 
we will report the simultaneous prediction of electron, 
spin, and lattice edge modes, which should further facili- 
tate the identification of the intrinsically inhomogeneous 
phases by cross correlating electronic gap modes that oc- 
cur on eV scale with the localized lattice and spin modes 
that appear on meV scales. In addition to ARPES and 
inelastic magnetic and lattice neutron scattering, exper- 
imental techniques that could be applied to study the 
edge modes include infra-red and Raman spectroscopy 
[24]; resonant Raman scattering would be particularly 
valuable even for low doping levels [25]. 
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